from sentinelsat import SentinelAPI, read_geojson, geojson_to_wkt


# 创建SentinelAPI，请使用哥白尼数据开放获取中心自己的用户名及密码
api =SentinelAPI('liangjiahao', 'liangjiahao', 'https://scihub.copernicus.eu/dhus/')

# 读入上海市的geojson文件并转换为wkt格式的文件对象，相当于足迹
footprint =geojson_to_wkt(read_geojson('Henan_map.geojson'))

# 通过设置OpenSearch API查询参数筛选符合条件的所有Sentinel-2L2A级数据
products =api.query(
	footprint,                        # Area范围
	date=('20200903', '20200910'),   # 搜索的日期范围
	platformname='Sentinel-2',        # 卫星平台名，Sentinel-2
	producttype='S2MSI2A',            # 产品数据等级，‘S2MSI2A’表示S2-L2A级产品
	cloudcoverpercentage=(0, 30))   # 云量百分比
		
# 通过for循环遍历并打印、下载出搜索到的产品文件名
for product in products:
	# 通过OData API获取单一产品数据的主要元数据信息
	product_info = api.get_product_odata(product)
	# 打印下载的产品数据文件名
	print(product_info['title'])
	# 下载产品id为product的产品数据
	api.download(product)
